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Abstract 

A forward amplitude analysis on pp and pp elastic scattering above 5 GeV 
is presented. The dataset includes the five recent high-precision TOTEM 
measurements of the pp total cross section at 7 and 8 TeV. Following previ- 
ous works, the leading high-energy contribution for the total cross section is 
parametrized as ln 7 (s/s^), where s is the cm. energy squared, 7 and Sh are 
free real fit parameters. Using singly-subtracted derivative dispersion rela- 
tions the total cross section (a to t) and the rho parameter (p) are analytically 
connected. Different fit procedures are considered, including individual fits to 
a to t data, global fits to cr tot and p data, constrained and unconstrained data 
reductions. The results favor a rise of the total cross section at the LHC 
energy region faster than the log-squared bound by Froissart and Martin, 
namely 7 greater than 2. A critical discussion on the correlation, practical 
role and physical implications of the parameters 7 and Sh, including a physi- 
cally meaningful representation for Sh, is presented. The discussion comprises 
the 2002 prediction by the COMPETE Collaboration and the recent result 
by the Particle Data Group (2012 edition of the Review of Particle Physics). 
Some conjectures on possible implications of a fast rise of the proton total 
cross section at the highest energies are also presented. 
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1. Introduction 

The theoretical and phenomenological descriptions of the energy depen- 
dence of the hadronic total cross-section at high energies have been a funda- 
mental and long-standing problem. Once the optical theorem connects the 
total cross section with the imaginary part of the forward elastic scattering 
amplitude, the main theoretical difficulty concerns the lack of a pure (model- 
independent) nonperturbative QCD description of the soft scattering states, 
in particular the elastic channel in the forward direction. In this introduc- 
tion, we shortly recall some results on the rise of the total cross section of 
interest in this work, followed by the plan of the paper. 
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Historical formal analyses by Froissart, Lukaszuk and Martin [l], 0, 0, 0] 
have established the famous asymptotic bound for the rise of the total cross 
section, 

0tot(s) < cln 2 — , 
so 

where s is the cm energy squared, c < vr/m 2 60 mb and So is a constant. 
After the Martin derivation in the context of the axiomatic quantum field 
theory jijj], the log-squared bound has played a central role in model construc- 
tions, aimed to treat, interpret and describe strong interactions. Recently, 
Ya. I Azimov has argued that it is not obvious if QCD can be considered 
an axiomatic theory and has also established new constraints for the rate of 
increase of crtot(s) pi- The formalism allows to conclude that, depending on 
the behavior of the amplitude in the nonphysical region, the total hadronic 
cross section may rise faster than the log-squared bound, without violation 
of unitarity (see also @, 0|). 

The dependence of the total cross section on the energy has been usually 
investigated through forward amplitude analysis. This approach is charac- 
terized by different analytical parametrizations for o~tot( s ) an d global data 
reductions including the p parameter (ratio between the real and imaginary 
parts of the forward amplitude), using singly-subtracted dispersion relations. 
In 2002 the COMPETE Collaboration has developed a detailed analysis on 
different functional forms within a ranking procedure, including several re- 
actions and different energy cutoffs js), [if. The analysis has favored the 
asymptotic form ln 2 (s/so), when compared with ln(s/s ) or powers [s/so] e 



with e above 1, a conclusion corroborated by other analysis [10J, lUj. After 
that, data reductions with the COMPETE highest-rank parametrization (in 
agreement with the log-squared bound and all the experimental data avail- 
able), became standard references in the Review of Particle Physics by the 



Particle Data Group (PDG) up to the 2010 edition (12 . 

In the experimental context, the recent results from the LHC, reaching 
the highest energy region, opened new expectations in the phenomenological 
and theoretical contexts. However, the first results on pp elastic, diffrac- 
tive and non-diffractive processes by the TOTEM Collaboration, have indi- 
cated noticeable discrepancies with standard predictions from representative 



phenomenological models [13l . Il4j . Il5| . In general grounds, this novel exper- 



imental information seems not to be yet able to select or exclude classes 
of phenomenological approaches and/or theoretical pictures in a clear and 
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conclusive way. 

In what concerns the pp total cross-section, the first 7 TeV TOTEM 
high-precision measurement has been obtained with small bunches and a 
luminosity-dependent method, indicating 

o~tot,i(7 TeV) = 98.3 ± 2.8 mb (small bunches / luminosity-dependent) flij ]. 

In particular and for our purposes, the COMPETE prediction (published 
10 years ago) is in good agreement with this datum [lj]. However, contrasting 
with this striking result, the 2012 edition of the Review of Particle Physics 
by the PDG has shown that with an updated dataset, including the above 
TOTEM datum, the fit with the highest-rank parametrization selected by the 
COMPETE Collaboration, presents no agreement with the TOTEM result 



!fj. Both results (COMPETE 2002 and PDG 2012) are illustrated in Figure 



[T] and will be discussed along the paper. 

Recently, Fagundes, Menon and Silva 17J,|l8) have revisited a parametriza- 
tion for the total cross section introduced by Amaldi et al. in the seventies, 
characterized by the exponent in the high-energy leading logarithm contribu- 



tion as a free real fit parameter |19j. Based on different fit procedures, vari- 
ants and data ensembles from pp and pp forward scattering at \fs > 5 GeV, 
the authors have shown that, once the 7 TeV TOTEM datum is included in 
the dataset, several statistically consistent solutions are obtained with the ex- 
ponent greater than 2, indicating, therefore, a faster-than-squared-logarithm 



rise for the total cross section 17, 18 



More recently, the TOTEM Collaboration has obtained four new high- 
precision measurements for the total cross section, three data at 7 TeV and 



one datum at 8 TeV [20, 1211, |22j|. The results and experimental conditions 
can be summarized as follows: 

o~tot,2(7 TeV) = 98.6 ± 2.2 mb (large bunches / luminosity-dependent) j2oj, 

(JtotA 7 TeV ) = "- 1 ± 43 mb (P independent) 0, 

o~tot,4:(7 TeV) = 98.0 ± 2.5 mb (luminosity-independent) (2l| . 

<Jtot(8 TeV) = 101.7 ± 2.9 mb (luminosity-independent) [22]. 

In this work we revisit, once more, the Amaldi et al. parametrization, 
now taking into account the above recent measurements of the pp total cross 
section at 7 and 8 TeV. The general strategy follows the approach developed 



in our previous analysis [17J, [18|. However, besides the updated dataset, the 



focus here is on the practical role and the physical meaning of the free fit 



4 



parameters associated with the high-energy leading contribution. In partic- 
ular, we discuss the differences between the 2002 COMPETE prediction and 
the 2012 PDG result in the phenomenological context (Reggeon/Pomeron 



exchanges). As in |17I . |l8j, the analysis is also restricted to pp and pp data 



at y/s > 5 GeV (the same energy cutoff used by the COMPETE Collabora 



tion |8|, |9[ and in the last PDG versions [12|, [16|). Our analysis and results 
favor, once more, a rise of the hadronic total cross section faster than the 
log-squared bound by Froissart and Martin, at the LHC energy region. 

The manuscript is organized as follows. In section [2] we recall the ana- 
lytical parametrization for o~ to t(s) and the formula connecting this quantity 
with p(s), using singly-subtracted derivative dispersion relations 0; the fit 
procedures and strategies are also outlined. In section [3] we present the fit re- 
sults with datasets up to \fs max = 7 TeV and \fs max = 8 TeV, in the cases of 
unconstrained and constrained data reductions, individual fits to cr tot {s) data 
and global fits to a to t and p data. A discussion on all the obtained results, 
with focus on the high-energy leading contribution parameters, is presented 
in section HJ The conclusions and some final remarks are the contents of 
section [51 



2. Analytical Parametrization and Fit Procedures 



In this section, after introducing the two formulas to be employed in our 
data reductions (parametrization for o~ tot (s) and the analytical expression for 
p(s) from derivative dispersion relations), we outline some important points 
on the fit procedures and strategies. 

2.1. Analytical Parametrization and Dispersion Relation Result 

We consider the analytical parametrization for the pp and pp total cross 
section introduced by Amaldi et al. in the seventies 19] and also employed by 
the UA4/2 Collaboration in the nineties [23}. It consists of two components, 
associated with low- (LE) and high-energy {HE) contributions: 



o-tot{s) = a LE (s) + a HE (s) 



(1) 



The first term accounts for the decreasing of the total cross section and the 
differences between pp and pp scattering at low energies and in this work it 
is expressed by 



o-le(s) 







-bi 
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+ ra 2 
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(2) 
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where r = -1 (+1) for pp (pp) scattering, Si = 1 GeV 2 is fixed, while ax, bx, 
a 2 and b 2 are free fit parameters. 

The second term accounts for the rising of the cross section at higher 
energies and is given by 



phe(s) = a + (3 In 7 — , 

Sh 



(3) 



where a, (3, 7 and Sh are free real parameters. 

For further reference, we briefly recall that, in the context of the Regge- 
Gribov theory, the decreasing (Jle(s) contribution is associated with Reggeon 
exchanges: 61 and b 2 correspond to the intercept of the trajectories and ai, 
a 2 to the Reggeon strengths (residues). The <Jhe(s) term simulates the rise 
of the total cross section and is associated with the Pomeron exchange. For 
example, for 7 = 1, the constant plus Ins terms correspond to a double 
pole at J = 1 and for 7 = 2 a triple pole (expressed by In 2 s, In s and the 
constant terms) [§). For our purposes it is important to note that for 7 = 2, 
parametrization (1-3) has the same analytical structure of the highest-rank 
parametrization selected by the COMPETE Collaboration [§, ^ . 

The analytical connection with the p parameter is obtained using singly- 
subtracted derivative dispersion relations in the operator expansion form 
introduced by Kang and Nicolescu [2H (discussed also in 25|). In terms of 
the parametrization (1-3) for pp and pp scattering, the analytical results for 



p(s) read [171, |l8 



P(s) 






s 


m 





+ A In 7 " 1 ( — 

Sh 



+ B In 7 " 3 f — 1 + C In 7 " 5 f — 1 +ra 2 cot 



Sh 



where K is the subtraction constant, 



Sh 



7T n ^,1 

A=-M, B = - 



/3 7 [7-l][7-2], 



C 



2 

15 



7T 

2J 



/3 7 [7-l][7-2][7-3][ 7 -4], 
and, as before, r = -1 (+1) for pp (pp) scattering. 
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2.2. Fit Procedures and Strategies 

The fit procedures and methodology, as well as the role and applicability 
of the subtraction constant in the derivative dispersion relation approach, 



have been detailed discussed in [17] and [18] . Here we recall and also outline 
some important points of interest in this work, in special the introduction of 
a representation for the parameter Sh- 

2.2.1. Dataset 

The input dataset for fits concerns only accelerator data on a to t and p 
from pp and pp scattering, covering the energy region from 5 GeV up to 8 
TeV. The energy cutoff is the same employed in the COMPETE and PDG 
analyses jgl, [lij]. The data below 7 TeV have been collected from the PDG 
database, without any kind of data selection or sieve procedure. Estimations 
of the pp total cross section from cosmic-ray experiments will be displayed 
in the figures as illustrative results. The TOTEM estimation for p at 7 TeV 



2 II ] is also displayed as illustration. All the references on these data and 



estimations can be found in 18 . 



2.2.2. Nonlinearity and Feedback Values 

Once the nonlinearity of the fit demands a choice for the initial values 



(feedbacks) for all free parameters 26|], different choices have been tested 



and discussed in our previous analyses 17, 18[ (see also 27]). Here, to 
initialize our parametric set, we consider only the values of the fit results 
present in the updated 2012 PDG version, obtained with the highest-rank 
COMPETE parametrization j^, [l6j]. This choice can be considered rather 
conservative in the sense that 7 has been fixed to 2 in the data reductions. 
However, it should be noted that the PDG and COMPETE analyses include 
different collision processes (meson-p, baryon-p, among others) and also tests 
on universatility. Therefore our dataset, restricted to pp and pp scattering, 
correspond to only a subset of the data ensemble employed in the global 
analysis by the COMPETE and PDG. Despite of this, we understand that 
with this choice for the feedbacks, we initialize our parametrization with a 
statistically and physically meaningful input, contributing with the search 
for the best fit (see section 12.2.41 below). 

The values of the feedback parameters to be used in our fits are displayed 
in the third column of Table [H For further reference in section HI we dis- 
play in Table [1] the values of the parameters obtained with the highest-rank 
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parametrization (Eos. (1-3) with 7 = 2) in both the 2002 COMPETE anal- 
ysis (Table VIII in |j) and the 2012 PDG version (Table 46.2 in [16]), which 
is based on an updated dataset including the first 7 TeV TOTEM datum. 
In the last case, the values of the parameters a± and ci2 correspond to the 
normalization of Eq. (2), namely si = 1 GeV 2 fixed (which is different from 
the normalization adopted in [16|). The corresponding results for a to t(s) and 
uncertainty regions (evaluated through propagation from the errors in Table 
[1]) are shown in Figured], together with the experimental information (in this 
figure only the first 7 TeV TOTEM measurement is displayed, as in (l6|). 
We shall return to these results along the paper. 

Table 1: Fit results through parametrization (1-3) with 7 = 2 from the COMPETE 
and PDG analyses. The parameters ai, 02, a and /3 are in mb, Sh in GeV 2 , b\, 62 are 
dimcnsionless (s/ = 1 GeV 2 ). 





COMPETE 2002 [9] 


PDG 2012 [16] 




42.1 ± 1.3 


46.07 


± 0.76 




0.467 ± 0.015 


0.462 


± 0.002 


a 2 


32.19 ± 0.94 


34.02 


± 0.63 




0.5398 ± 0.0064 


0.550 


± 0.005 


a 


35.83 ± 0.40 


34.71 


± 0.15 





0.3152 ± 0.0095 


0.265 


± 0.050 


Sh 


34.0 ± 0.54 


16.21 


± 0.16 



2.2.3. Individual and Global Fits 

Global fits to o~tot and p data demand the use of dispersion relations with 
one subtraction and therefore the introduction of one more free parameter, 
the subtraction constant. As already discussed in 17j and jl8|] this parameter 
does not have a physical interpretation as those present in the total cross 
section parametrization (associated with Reggeon and Pomeron exchanges). 
Moreover, several authors have criticized the usual methods to extract p (see 
section 2.2 in [l7j]) and, in addition, due to the correlation among all the 



dobal fits constraints the rise of the total 
and references therein. Despite 



28, 29 



fit parameters, the p inclusion in 
cross section, as demonstrated in 
these disadvantages, here we shall consider both individual fits to a to t data 
through Eqs. (1-3) and global fits to a to t and p data using Eqs. (1-5). 
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2.2.4- Minimization and Statistics 

The data reductions have been performed with the objects of the class 
TMinuit of ROOT Framework (Sj. We have employed the default MINUIT 



error analysis [3JJ with the selective criteria that follow. In the minimization 
program a Confidence Level of one standard deviation was adopted in all fits 
(UP = 1). In each test of fit, successive running of the Migrad have been 
considered (up to 5,000 calls), until full convergence has been reached, with 
the smallest FCN (chi-squared) and EDM (Estimated Distance to Minimum), 
specifically EDM < 10~ 4 (adequate for the one sigma CL). In addition, the 
error in the parameters should not exceed the central value. The error matrix 
provides the variances and covariances associated with each free parameter, 
which is used in the analytic evaluation of the uncertainty regions in the fitted 



and predicted quantities, through standard error propagation procedures 26 

As test of goodness of fit we shall consider the chi-square per degree of 
freedom (x 2 /DOF) and the corresponding integrated probability, P(x 2 ) [26]. 
The goal is not to compare or select fit procedures or fit results but to check 
the statistical consistence of the data reductions in a reasonable way. 

2.2.5. Unconstrained and Constrained Fits 

To consider the exponent 7 in the leading logarithm component as a 
real (not integer) free fit parameter leads to some special consequences and 
particular conditions. These aspects, including the effects of the correlation 
between the free parameters 7 and Sh in data reductions, will be discussed in 
some detail in section HJ after the presentation of our fit results. In order to 
implement and facilitate that discussion, we introduce here a useful represen- 
tation for the high-energy scaling factor Sh, which will led us to distinguish 
between unconstrained and constrained fits, as explained in what follows. 

The representation is based on two arguments: (i) the reference to Sh as 
a (unknown) constant in the Froissart-Martin derivation of the bound; (ii) 
the reasonable physical conjecture that this factor might be proportional to 
the energy threshold for the scattering states (above the resonance region), 
namely Sh oc (m p + m p ) 2 , where m p is the proton mass. In this case, we can 
represent the scaling factor by 

s h = 5 [4m 2 p ], Am 2 p = 3 ' 521 Gey2 ' ( 6 ) 

where 5 is a real dimensionless parameter and 5 = 1 at the physical threshold. 

With this representation, we can distinguish two physical conditions in 
data reductions: either to consider 5 indeed as a free fit parameter (equiv- 
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alently, Sh a free fit parameter) or to assume 5 = 1 (equivalently, to fix 
Sh = 4=iTip, the energy threshold). In what follows we will denote these two 
cases by unconstrained fits (5 free) and constrained fits (5 = 1 fixed). 

Here, as in the COMPETE and PDG analyses, we shall treat Sh as a 
constant (with the above representation). However, it should be noted that 
the possibility of a slow rise of Sh with s for large s has been discussed by 
some authors, as for example in jiij an also by Azimov Jsf] . 

2.2.6. Ensembles and Feedbacks 

As mentioned earlier, our dataset consists of pp and pp accelerator data 
at > 5 GeV. In order to investigate the effect in the fits associated with 
the TOTEM measurement at 8 TeV, as compared with the previous results 
at 7 TeV, we shall consider two data ensembles. The first one with data up 
to 7 TeV (including the four measurements) and a second one adding the 8 
TeV datum. For reference, we will denote these two cases by \fs max = 7 TeV 
ensemble and \fs max = 8 TeV ensemble, respectively. 

For the \fs max = 7 TeV ensemble, we use as feedback the values of the 
parameters from the 2012 PDG version, displayed in the third column of 
Table CD In the case of global fits to a to t and p data, we consider K = for 
the initial value of the subtraction constant [16]. The fit results with this 
ensemble are then used as feedback to initialize the parametrization with the 
\fs max = 8 TeV ensemble, in each one of the variants considered (individual, 
global, unconstrained and constrained). 

3. Fit Results 

Summarizing the discussion in the last section, we select two ensembles 
of accelerator data on a tot and p at y/s > 5 GeV, one up to \fs max = 7 
TeV and another one up to \fs max = 8 TeV. In each case we consider both 
unconstrained fits (st or 5 in Eq.(6) as free fit parameter) and constrained 
fits (sh = ^rn 2 or 5 = 1 fixed), treating also both individual fits to a to t data 
through Eqs. (1-3) and global fits to a to t and p data, using Eqs. (1-5). In 
what follows we present the fit results and discuss the accordance with the 
selective criteria outlined in section 12.2.41 

3.1. yfs max = 7 TeV Ensemble 

For this ensemble all data reductions presented agreement with the se- 
lective criteria. The fit results and statistical information are displayed in 
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Table [2J The curves, uncertainty regions (from error propagation) and ex- 
perimental information in the case of global fits to a to t and p are shown in 
Figure [2] for the unconstrained fit and in Figure [3] for the constrained one. 



Table 2: \fs max = 7 TeV Ensemble. Fit results through Equations (1-3) for the a to t data 
and Equations (1-5) for the atot and p data in the case of unconstrained (s/, free) and 
constrained (sh = i-nri) data reductions. Units as in Tabled] 





Unconstrained (5 free) 


Constrained (5 = 1) 




O~tot 


a to t and p 


O~tot 


a to t and p 


ax 


59.5 ± 8.8 


56.7 ± 6.4 


60.5 ± 7.8 


59.1 ± 4.9 


h 


0.553 ± 0.063 


0.528 ± 0.062 


0.526 ± 0.074 


0.500 ± 0.054 


a 2 


33.2 ± 2.2 


34.0 ± 2.0 


33.2 ± 2.3 


34.1 ± 1.9 


b 2 


0.541 ± 0.015 


0.547 ± 0.013 


0.541 ± 0.016 


0.547 ± 0.013 


a 


35.53± 0.70 


35.3± 1.0 


34.0 ± 1.6 


33.2 ± 1.3 


(3 


0.184 ± 0.066 


0.242 ± 0.064 


0.106 ± 0.044 


0.130 ± 0.042 


7 


2.14 ± 0.12 


2.053 ± 0.094 


2.29 ± 0.16 


2.21 ± 0.11 


Sh 


12.8 ± 3.0 


16.41 ± 0.64 


3.521 (fixed) 


3.521 (fixed) 


K 




40 ± 14 




39 ± 12 


DOF 


159 


234 


160 


235 


X 2 /DOF 

p(x 2 ) 


0.92 


1.09 


0.91 


1.09 


0.754 


0.155 


0.775 


0.164 


Figure: 


2 


3 



3.2. ^s max = 8 TeV Ensemble 

For this ensemble, the unconstrained fits presented some disagreement 
with the selective criteria. In the case of global data reduction to a to t and p 
the fit did not reach the full converge (the fit failed) and therefore we have 
no solution for this case. In the individual fit to a to t the data reduction 
converged but it should be noted that the corresponding error matrix is not 
positive definite. 

The fit results and statistical information are displayed in Table El The 
curves, uncertainty regions and experimental information in the case of the 
constrained global fit to atot and p are shown in Figure H] (the unconstrained 
case did not converge). 
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Table 3: \fs max = 8 TeV Ensemble. Fit results through Equations (1-3) for the a to t data 
(unconstrained, Sh free, and constrained, Sh = 4nip) and Equations (1-5) for the er tot and 
p data (only the constrained case) . Units as in Table [T] 





Unconstrained (S free) 


Constrained (5=1) 




o'tot vtot and p 




a to t and p 




61.4 ± 1.7 


60.8 ± 6.9 


59.5 ± 7.2 


h 


0.5062 ± 0.0084 


0.530 ± 0.061 


0.505 ± 0.075 


a 2 


33.2 ± 1.7 


33.2 ± 2.3 


34.1 ± 2.0 


b 2 


0.541 ± 0.012 


0.541 ± 0.016 


0.547 ± 0.013 


a 


32.37 ± 0.22 


34.1 ± 1.2 


33.4 ± 1.8 





0.0477 ± 0.0022 


0.102 ± 0.033 


0.124 ± 0.054 


7 


2.493 ± 0.016 


2.30 ± 0.11 


2.23 ± 0.15 


Sh 


0.633 ± 0.073 


3.521 (fixed) 


3.521 (fixed) 


K 






40 ± 16 


DOF 


160 


161 


236 


x 2 /dof 
p(x 2 ) 


0.91 


0.91 


1.09 


0.774 


0.788 


0.172 


Figure: 


5 (up) 


5 (down) 


4 



4. Discussion 

Using parametrizations (1-5), we are interested in a consistent quantita- 
tive description of the rise of the total cross-section at high energies. We 
have considered 8 variants of data reductions, obtaining full convergence in 
7 cases. The fit results are displayed in Tables |2] and [3J Our main goal is 
to investigate if these results indicate a log-squared behavior or a rise faster 
than this bound. In what follows, for a given result on 7± A7, we consider a 
result statistically consistent with a faster rise the cases in which 7 — A7 > 
2. In this case, for short, we will refer to a result statistically consistent with 
7 > 2. However, as remarked early in section I2.2.5[ to consider the expo- 
nent 7 as a real (not integer) fit parameter has some special implications not 
present in the canonical assumption 7 = 2 (fixed). These aspects are directly 
related to the high-energy scaling factor Sh and have important consequences 
not only on the data reductions but also on the physical interpretation of the 
fit results. To discuss these aspects is the goal of this section. 

To this end, we first present our partial conclusions on the variants and 
fit results, with focus on the value of the parameter 7 and its relation with 
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a faster-than-squared-logarithm rise. After that we discuss the practical role 
and physical implications associated with the correlation between 7 and s^. 
This discussion will led us to a final conclusion in favor of the constrained 
fits as our best results, in both physical and statistical contexts, indicating 
a rise of the total cross section faster than the log-squared bound. At last, 
we present some conjectures on possible implications of this conclusion. 

4-1. Partial Conclusions on the Fit Results 

Based on the results displayed in Tables [2] and [3j we present a compar- 
ative discussion related to two aspects: individual/global fits and uncon- 
strained/constrained fits. 

4-1.1. Individual and Global Fits 

From Tables M(y/s max = 7 TeV Ensemble) and[3](y / s maa . = 8 TeV Ensem- 
ble) all the individual fit results to a tot data (constrained and unconstrained) 
are statistically consistent with 7 > 2 (confirming, therefore, the conclusions 
first presented in [I3|). In all cases the integrated probability reads P{x 2 ) ~ 
0.8. The highest 7-values are associated with the \fs max = 8 TeV Ensemble 
(the largest energy interval with experimental information presently avail- 
able), indicating 7 rs 2.5 (unconstrained fit) and 7 ~ 2.3 (constrained fit). 
The corresponding results for a to t(s) and uncertainty regions are displayed 
in Figure |5] for both the unconstrained and constrained fits, together with 
the experimental data. We notice that the agreement with the TOTEM 
measurements at 7 TeV is striking. In the case of the constrained fit the 
uncertainty region reaches the central value at 8 TeV including also all the 
four 7 TeV central values. The numerical results and predictions for the total 
cross section at some energies of interest are displayed in Table H] from all the 
variants investigated in the individual fits. We also note that with the \fs max 
= 8 TeV Ensemble, the predictions at 57 TeV, namely atot ~ 142 - 143 mb 
are about 7 % larger than the central value of the result by the Pierre Auger 



Collaboration, namely 133 mb [33 



In the case of global fits to a tot and p data the results depend on the 
ensemble and on the constraint condition considered. In all cases of global 
convergent fits the integrated probability reads P(x 2 ) ~ 0.2. For the con- 
strained fits (both ensembles) the results are statistically consistent with 7 > 
2. In the unconstrained case and \fs max = 7 TeV Ensemble the result may 
be considered barely consistent with 7 > 2, since, up to 2 figures, 7 lies in 
the interval 2.0 - 2.2. For the \fs max = 8 TeV Ensemble we did not obtain 



13 



Table 4: Fit results and predictions for the pp total cross section at higher energies from 
individual fits to a to t data. 





v^ max = 7 TeV 
Unconstrained Constrained 


Ana* = 8 TeV 

Unconstrained Constrained 


7 TeV 

8 TeV 
14 TeV 
57 TeV 


97.8 ± 1.1 97.9 ± 1.3 
100.2 ± 1.2 100.2 ± 1.4 
110.6 ± 1.6 110.8 ± 2.0 
140.6 ± 3.2 141.3 ± 4.4 


98.12 ± 0.86 98.1 ± 1.2 
100.56 ± 0.91 100.5 ± 1.3 
111.4 ±1.1 111.2 ±1.7 
142.7 ± 1.8 142.0 ± 3.5 



full convergence. The numerical results and predictions for the total cross 
section at the energies of interest are displayed in Table [5j We note that at 
57 TeV the results indicate a to t ~ 139 - 140 mb, which is about 5 % larger 
than the Auger central value. 



Table 5: Fit results and predictions for the pp total cross section at higher energies from 
global fits to a tot and p data. 





Vs m ax = 7 TeV 
Unconstrained Constrained 


v^ = 8TeV 
Unconstrained Constrained 


7 TeV 

8 TeV 
14 TeV 
57 TeV 


97.5 ± 1.2 97.6 ± 1.2 
99.8 ± 1.3 99.9 ± 1.3 
109.9 ± 1.7 110.2 ± 1.7 
138.8 ± 3.3 139.6 ± 3.4 


97.8 ± 1.3 
100.2 ± 1.4 
110.6 ± 2.0 
140.4 ± 4.3 



From the above discussion on the 7 value, we conclude that our results fa- 
vor an increase of the hadronic total cross-section faster than the log-squared 
bound at the LHC energy region. The values of the 7 parameter from all the 
full convergent fits, within the uncertainties, are schematically displayed in 
Figure [6j 

4-1-2. Unconstrained and Constrained Fits 

In this respect we note a substantial negative correlation between the 
parameters 7 and Sh- With the \fs max = 7 TeV ensemble, in going from the 
unconstrained to the constrained case, 7 increases and Sh decreases (in both 
individual and global fits; with the \fs max = 8 TeV ensemble (individual fit) 
7 decreases and Sh increases. These effects are associated with the correlation 
among all the fit parameters due to the nonlinear data reductions. 
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In the MINUIT code the symmetric correlation matrix provides a measure 
of the correlation between each pair of free parameter through a coefficient 



with numerical limits ± 1 (full correlation) and (no correlation) [26, 30| . For 
example, the coefficients in the global fits with the \fs max = 7 TeV ensemble, 
in the cases of the unconstrained and constrained fits, are displayed in Table 
El 

We also note from Figure El that in the case of constrained fits the values 
of the 7 parameter, 2.2 - 2.3, are more stable than in the unconstrained 
cases. We understand that all these characteristics favor the constrained fits. 



Table 6: Correlation coefficients from the correlation matrices associated with uncon- 
strained fit (UF) and constrained fit (CF) in the case of global fits with the \fs max = 7 
TeV ensemble. The off-diagonal coefficients from the UF are displayed above the diagonal 
of the table (not filled) and those from the CF, below that diagonal. 







UF 






ai 


h 




b 2 


a 


P 


7 


Sh 


K 




a\ 




0.981 


0.154 


0.149 


0.933 


-0.904 


0.878 


-0.598 


0.914 




bi 


0.966 




0.092 


0.094 


0.985 


-0.969 


0.950 


-0.539 


0.909 






0.193 


0.129 




0.985 


0.058 


-0.052 


0.054 


0.183 


-0.018 


CF 


b 2 


0.178 


0.122 


0.985 




0.063 


-0.059 


0.062 


0.187 


-0.038 


a 


0.909 


0.985 


0.105 


0.100 




-0.997 


0.988 


-0.443 


0.877 




P 


-0.847 


-0.949 


-0.125 


-0.122 


-0.986 




-0.996 


0.389 


-0.855 




7 


0.825 


0.933 


0.126 


0.123 


0.976 


-0.998 




-0.320 


0.833 




Sh 


















-0.543 




K 


0.882 


0.872 


0.000 


-0.032 


0.835 


-0.787 


0.770 







4-2. The Role and Effects of the Parameters 7 and Sh 

Based on the results and discussion in the last section, we treat here the 
important practical and physical role of the scale parameter Sh, specially in 
data reductions with 7 as a free real (not integer) parameter. We first list five 
characteristics of interest involved, distinguishing the case of 7 = 2, and then 
discuss the connections with our fit results, as well as with the COMPETE 
2002 and PDG 2012 results (Table [Hand Figured]). 

1. Even treating Sh as an unknown constant, with the canonical choice 7 
= 2 the scaling factor can be estimated (or determined) through the data 
reduction. However, in the case of 7 real and free we have two unknown and 
correlated parameters 7 and s/,. As a consequence, we can obtain solutions 
statistically consistent for different values of these parameters and that may 
imply in different physical pictures. 
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2. For 7 = 2 (fixed) the high-energy contribution (The( s ), Eq. (3), is well 
defined for all values of s as compared with the Sh value. The only difference 
concerns the fact that for s > Sh the logarithm term increases with the energy 
and for s < Sh this contribution decreases as the energy increases (we will 
return to this important point in what follows). On the other hand, for 7 
real (not integer) the logarithmic term is not defined for s < Sh so that this 
component can only starts at s = Sh- Above this point the contribution 
increases with the energy. 

3. For 7 = 2 and Sh = 8 [4m^] the high-energy contribution can be expressed 

by 



,„(.) = . +11b (_J + ^_|, (7) 

where a = a + /31n 2 (<5) and b = — 2/31n(<5), providing an explicit correlation 
among the high-energy parameters. The above expansion, however, is not 
possible in the case of 7 real (7^ 2) so that the correlations are somewhat 
hidden in the nonlinear data reduction. 

4. As commented before, the <Jhe{s) component accounts for the rise of 
the total cross section at high energies, being associated with the Pomeron 
exchange. In the standard soft Pomeron concept this term is expected to 
increase with the energy, as is the case of the simple pole parametrization, 



s e with e slightly greater than 1 3J, |35|, |36|, [37 



5. Another aspect of interest concerns the minimum value of the energy 
above which a given parametrization is supposed to be applied, or valid. In 
this respect, as we will discuss, the energy cutoff for data reduction, yfs min 
plays a central role, in connection with the corresponding energy scale, ^/sh. 
Here, as in both COMPETE and PDG analyses, we have adopted \fs min = 
5 GeV. 

Based on the above characteristics, it is expected that, depending on the 
values of s min and Sh, different physical interpretations can emerge in the 
cases of 7 = 2 (fixed) and 7 real (not integer). Let us discuss these aspects 
and their connection with the fit results. 

• 7 = 2 (fixed) 

In this case, if Sh < s m i n then, in the region of experimental data (y/s > 
y/s min ), the high-energy component <The(s) increases with the energy, as 
expected in the standard concept of the soft Pomeron contribution. That is 



16 



the case of the 2012 PDG version since, from Table [TJ Sh = 16.21 ± 0.16 < 
s mm = 25 GeV 2 . 

On the other hand, if Sh > s m i n then, in the interval \fs min < \fs < 
y/~Sh, the cthe(s) component decreases as the energy increases, suggesting a 
disagreement with the standard Pomeron contribution. That, however, is 
the case of the 2002 COMPETE result, since, from Table H s h = 34.00 ± 
0.54 > s min = 25 GeV 2 . The dependence of cthe(s) for the above two cases 
is illustrated in Figure [7] in the energy interval 4 GeV < \/s < 7 GeV that 
includes the energy cutoff \/s-„ = 5 GeV. 

t ~>*j v i lit 1 1 I I 

In this respect, according to the COMPETE Collaboration [8] : "One must 
note that in some processes, the falling ln 2 (s/so) term from the triple pole at 
s < So is important in restoring the degeneracy of the lower trajectories at low 
energy. Hence the squared logarithm manifests itself not only at very high 
energies, but also at energies below its zero." However, even accepting this 
argument on a decreasing Pomeron contribution, the fast divergent rise of 
&he{s) as the energy decreases below ^/sh (Figure [7]) and above the physical 
threshold, remains, in our opinion, unexplained. 

A practical or even pragmatic consequence of the COMPETE result is 
that the agreement between their 2002 prediction (with 7 = 2 in accordance 
with the Froissart- Martin bound) and the 7 TeV TOTEM measurement, is 
directly connected with the rather large value of s/,. By contrast, the PDG 
result with 7 = 2, smaller Sh and consistent with the rise of <Jhe{s) in the 
whole interval of energy investigated, lies below the TOTEM datum (compare 
Figures [1] and [7j) . 

Summarizing, with 7 = 2: (a) the COMPETE correctly describes the 7 
TeV TOTEM datum, but with a decreasing (The{.s) contribution above the 
cutoff s min up to Sh and a divergent (increasing) contribution as the energy 
decreases below s min ; (2) in the PDG 2012 result, cr HE (s) increases in the 
whole interval of energy investigated (above the cutoff) but lies (within the 
uncertatiny) below the TOTEM datum 

• 7 real (not integer) 

For 7 not integer, as commented before, the <The(s) component is not 
defined at \/s < s/s~h. The contribution starts at y/s = with ghe{.S}i) 
= a and from this point on it increases with the energy, as expected in the 
standard soft Pomeron exchange. In respect to our fit results, the <the(s) 
component is well defined in the whole interval of energy investigated since 
\fs min = 5 GeV and in all data reductions Sh < s min (Tables |2] and [3]) . 
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At last, based on the above physical aspects related to the scale factor, 
we discuss our unconstrained and constrained fit results in connection with 
our representation (6). This discussion will also led us to conclude that the 
best results are those obtained with the constrained case. 

With the unconstrained fits (5 or Sh = 54m 2 free) we have obtained for 
the \fs max = 7 TeV ensemble Sh ~ 13 GeV 2 (individual fit), Sh ~ 16 GeV 2 
(global fit) and for the \/s max 8 TeV ensemble Sh ~ 0.6 GeV 2 (individual 
fit). In these cases it seems difficult to devise a physical meaning for the 
beginning of the o~he{.s) component, because its value depend on the data 
analyzed (Figure |6]). Moreover, as mentioned in section |3T2| with the \/s max 
8 TeV ensemble, the error matrix is not positive definite in the individual fit 
and no convergence was obtained in the global fit. 

On the other hand, for the constrained fits (sh = 4m 2 fixed) the cthe(s) 
component starts at this threshold with aHE^Sh) = a and from this point on 
it increases with the energy, namely cthe(s) > a at yfs > 2m p . This situation 
seems to us physically meaningful in both phenomenological and theoretical 
contexts. Furthermore, in all cases investigated with the constraint condi- 
tion we have obtained full convergence in the data reductions and consistent 
statistical results. The values of 7 are also consistent (stable), around 2.2 - 
2.3, in all cases investigated. Therefore, we conclude that the constrained fits 
represent our best results, indicating a rise of the total cross section faster 
than the log-squared. 

4-3. Some Conjectures on a Fast Rise of the Total Cross Section 

At this point, we could conjecture (if not speculating) on the implication 
of a possible increase of a tot faster than In 2 s. One possibility points to a 
power-like behavior s e , e > 1 , which has always been an important and 



representative approach [38|, |39|, |4CJ, |4l|] . Predictions from unitarized models, 



developed near 10 years ago and consistent with the first 7 TeV TOTEM 



datum, are discussed, for example, in [42j, |43|, |44|, |45 



A faster-than-squared-logarithm rise could also indicate the onset of some 
new physics effect at the LHC energy region. For example, one possible ex- 
planation for the short penetration depth, recently observed in ultra-high- 
energy cosmic rays (UHECRs) around 100 TeV, is just an increase of the 
proton cross section faster than the extrapolations from models, which have 



been tested only at lower energies (see, for example, |46l . 1471 ] and references 
therein). These conjectures are not in disagreement with the recent theoret- 
ical arguments by Azimov 0,0, a. 
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At last, in contrast with an effective violation of the Froissart-Martin 
bound, a fast rise of the total cross section may also be associated with some 
local effect at the LHC energy region and/or beyond, so that, asymptotically, 
the bound might remain valid. In that case, 7 could represent a kind of 
effective exponente, depending on the energy and, possibly, associated with 
sums of different high-energy contributions. 



5. Conclusions and Final Remarks 

In 2002 Barone and Predazzi stated (High- Energy Particle Diffraction, 
page 140) 0: 

"The issue of the exact growth with energy of the total cross 
sections is both delicate and unresolved; the mild growth of total 
cross sections could be simulated by essentially any form and 
logarithmic physics is exceedingly difficult to resolve in a clear 
cut way." 



The aim of this work has been to take one more step in our investigation 
on the rise of the hadronic total cross-section at high energies. As in our 
previous analyses 17J, |l8j we have employed the analytical parametrization 
introduced by Amaldi et ai, with the exponent 7 in the high-energy leading 
logarithm contribution treated as a free real parameter in nonlinear data 
reductions. 

Here, the main points consisted in an updated analysis (including in the 
dataset the recent high-precision TOTEM measurements at 7 and 8 TeV) and 
a discussion on the correlation, practical role and physical meaning associ- 
ated with the exponent 7 and the energy scale factor Sh- As in our previous 
works, we have considered different variants, involving two ensembles, indi- 
vidual/global fits and unconstrained/constrained fits. As feedbacks for the 
nonlinear data reductions we have used the "conservative" values obtained 
by the PDG in the recent 2012 Review of Particle Physics edition. 

In section H] we have discussed all the fit results, indicating the advantages 
of the constrained fits (sh = 4m^ fixed) in both phenomenological and theo- 
retical contexts and noticing also the statistical consistence of the fit results. 
In special, we have concluded that the constrained fits with the \fs max = 8 
Te V ensemble (individual and global cases) represent our best results (Table 
[31 fourth and fifth columns and Figures 4 and 5). They indicate a rise of the 
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hadronic total cross-section faster than the log-squared bound at the LHC 
energy region. The results and predictions for the pp total cross section at 
energies of interest are displayed in the last columns of Tables H] and |5j 

We have conjectured that this fast rise might be connected with the onset 
of some new phenomena and have also speculated on the possible connection 
with the short penetration depth recently observed in UHECRs. If these 
effects have a local character (finite energies), there seems to be no contra- 
diction with the Froissart-Martin bound, since it has been derived for the 
asymptotic energy limit, s — > oo. 

In what concerns our approach, two aspects seems important to be ad- 
dressed. One of them consists on the extension of the parametrization to 
other physical quantities and the investigation of the physical consequences, 



as done in 18]. Another one concerns the possibility that the energy scale 
Sh might itself slowly grow with the energy j^, [32[. We are presently inves- 
tigatingthe above two aspects. Our final conclusion, as stressed in previous 



works [17|, [18j , is that the rise of the hadronic total cross-section at the high- 
est energies still constitute an open problem, demanding, therefore, further 
investigation. 
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Figure 1: Results with the COMPETE highest-rank parametrization (Eqs. (1-3) with 7 = 
2) from the 2002 COMPETE analysis and from the 2012 PDG version [l6| (which 

includes the first 7 TeV TOTEM datum) . The corresponding values of the parameters are 
displayed in Table Q] 
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Figure 2: Global unconstrained fit results with the \fs max — 7 TeV ensemble (third column 
in Tabled]). 



25 



• accelerator pp 
o accelerator pp 
t ARGO-YBJ 

a Akeno 

■ Fly's Eye 

* Auger 

♦ TOTEM 



[ J 











r 



_l I III' 



10 



10* 



10 J 
(GeV) 



10 4 



10° 




Figure 3: Global constrained fit results with the \fs max = 7 TeV ensemble (fifth column 
in Tabled 
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Figure 4: Global constrained fit results with the y/s max = 8 TeV ensemble (fifth column 
in Table EJ. 



27 



160 



• accelerator pp 
o accelerator pp 
t ARGO-YBJ 

a Akeno 

■ Fly's Eye 

* Auger 

♦ TOTEM 




unconstrained 



_l l iii' 



10 



10* 



10 J 
(GeV) 



10 4 



10= 



• accelerator pp 
o accelerator pp 
t ARGO-YBJ 

a Akeno 

■ Fly's Eye 

* Auger 

♦ TOTEM 




constrained 



_i I iii' 



10 



10* 



10 J 
(GeV) 



10 4 



10° 



Figure 5: Individual fit results to a to t data with the \/s max = 8 TeV ensemble and 
unconstrained (up) and constrained (down) data reductions (second column and fourth 
column in Tabled respectively). 
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Figure 6: Results here obtained for the exponent 7 as a free parameter in different data 
reductions through parametrization (1-3) and (1-5): y/s max — 7 and 8 TeV ensembles, 
constrained and unconstrained fits, individual (a to t) an d global (a to t and p) fits. 
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